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Abstract 

A quadratic invariant is denned as a quadratic form in the elements of a tensor that remains 
invariant under a group of coordinate transformations. It is proved that there are 7 quadratic 
invariants of the 21-element elastic m odulus tenso r un der SO (3 ) and 35 under SO(2). This 
answers some open questions raised by Ting (1987) and Ahmadl ( 2002 ) in this journal. 

-i-S ■ 

£ 1 Introduction 

-a . 

Ch The tensor of elastic moduli CijM is known to possess two linear invariants under arbitrary proper 
O orthogonal coordinate transformations, or SO (3): 

— Cijij — Cn + C22 + C33 + 2(C44 + C55 + Cqq), 
A = (kijj = Cn + C 22 + C33 + 2(ci2 + C 23 + C13). (1) 

1^ lAhmad ( 2002f) presented four quadratic invariants under SO(3), the first two of which were reported 

(NbylllnauS, 

7— 1 ■ 

^ ■ B\ = CijklCijkl, B 2 = CuklCjjkU -B3 = CukiCjkjU -B4 — CkiilCkjjl- (2) 

"£3 Ahmad demonstrated that the seven quadratic invariants {Af, A1A2, B±, B2, B3, B^} are in- 
£h dependent but did not show completeness. It is clear that the following is an eighth invariant, 

B5 = CijkiCikju (3) 

q although it is not so obvious whether or not it is independent of the other seven. We will prove 
that there are at most seven independent quadratic invariants under SO (3), and that the seven 
•'-j identified by Ahmad form a complete basis. In particular, B§ is a linear combination of this basis, 
rN specifically (see the Appendix) 

B s = -A\ + -Al - A,A 2 + B, - 2B 2 + 4B 3 - 2B t . (4) 
Consequently, every fourth order elasticity tensor satisfies this identity: 

It is well known that there are 5 linear invariants under rotation about an axis, or SO (2). Taking 
the axis as e 3 these are 

L\ = Cn + C22 + 2c 66 , L 2 = C44 + C55, L 3 = Cn + c 2 2 + 2ci2, £4 = C13 + C23, L 5 = C33. (6) 
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Ahmadl (120021 ) listed 17 quadratic invariants for SO (2): 



E2 
E 3 
£4 
E 5 
E 6 
E 7 
E$ 
E 9 
E10 
E u 

E\2 

En 
Eu 

E15 

Em 
En 



c 34 + c 35' 

(C15 + C 25 ) 2 + (C14 + C 24 ) 2 , 

(C15 + C 46 ) 2 + (c 24 + C 56 ) 2 , 

(C14 + C 24 )C 3 4 + (C15 + C 2 5)C 3 5, 

(C15 + 046)035 + (C24 + C5g)c 3 4, 

(C15 + C 25 )(ci 5 + C 46 ) + (C M + C 24 )(C 24 + C 56 ) 



"13 



~ l~ Coo ~~\~ *2iC 



-44 



,,2 
23 

,,2 
55 



+ 4 + 2c 



2 

36' 
2 

45' 



(cn + c i2 ) 2 + (c 12 + c 22 ) 2 + 2(ci fl + c 26 ) 2 , 
C13C55 + c 23 c 44 + 2c 36 c 4 5 

(en + ci 2 )ci 3 + (cia + c 22 )c 23 + 2(c 16 + c 26 )c 36 , 

(Cn + Ci 2 )c55 + (C12 + C 22 )C44 + 2{c\q + C 2 g)c45, 

4 + 4 + 24 + 4(4 + 4 + 4), 



(7) 



"14 



"I - C 2 4 + Cis + C 



15 



"25 



2(c 



06 



"46 V' 



(C44 - C 55 )c 36 + (C13 - C 23 )c 45 , 

(Cl3 — c 23)( c 16 + c 26) — ( c ll — C 22)C36, 

(C44 — C55)(Ci6 + C 2 6) + (Cxi — C 22 )C45- 



Tingl ( 1987 ) had previously reported 15 invariants, and these can be shown ( lAhmadl . 12002 ) to be 



contained in Ahmad's larger set. Ahmad also demonstrated that the 32 quadratic invariants formed 
from the 15 quadratic combinations of L%, . . . , L5 plus the 17 invariants E\, . . . , En are independent 
of one another. The question of completeness remained open. 

The purp ose of th i s pap er is to finish the study initiated by iTingl ( 19871 ) and subsequently 
expanded by lAhmadl (j2002l ). Two principa l resu lts are derived: first that the seven quadratic 
invariants under SO(3) identified by lAhmadl ( 20021 ) are indeed complete, and hence any quadratic 
isotropic invariant must be a linear combination of these. Secondly, we prove that there are 35 
quadratic invariants under SO(2). A compl ete basis for t he 35-dimensional space of quadratic 
invariants under SO (2) is formed by the 32 of I Ahma"dl (j2002l ) augmented by the these three, 

-^18 = ( c 15 + C46) c 14 — (c 2 4 + C^q)c25 — C15C56 + C 2 4C46, 

E19 = (C15 + C46)c 3 4 — (c 2 4 + 0^)035, (8) 
E20 = (Cl5 + C 25 )C 34 - (C14 + C 24 )C 35 . 

We consider quadratic forms on Ela, the space of tensors of elastic modulfl which possess the 
underlying symmetries 

Cijki — Cjiki, Cijki = Ckuj, (9) 
implying 21 independent elements, at most. The general quadratic form on Ela is 

EijklpqrsC*ijklCpq rsy (10) 

where Fijki pqrs are the elements of an eighth order tensor. Based on the properties (Q of the moduli, 
the elements of F satisfy 



ijklpqrs 



jiklpqrsi 



ijklpqrs 



F, 



klijpqrsi 



ijklpqrs 



pqrsijkl- 



lr ro be precise, Ela comprises tensors that are positive definite, but that distinction is not necessary here as we 
are concerned with quadratic forms on the moduli not defined by the moduli. 
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Identifying quadratic invariants is therefore equivalent to describing the properties of the eighth 
order tensor F, in particular the number of independent elements that survive under the group 
of transformations considered. The problem is linked to that of finding the integrity basis for 
fourth order tensors. As noted in a review on tensor functions: "Relatively little is known about 
representations of functions of tensors of order higher than two, for a ny of the transformation groups 
of interest in continuum mechanics" ( Rychlewski and Zhangi . Il99ll . p. 83). Considerable work has 
been done on this topic fo r secon d ord er tensors si nce this is critical to the form of elastic strain 
energy functions ( Spencerl . Il97ll ) (see IXiaol ( 19961 ) for recent developments on n=2). However, 
this literature is not directly applicable since we are concerned with tensors of higher order than 
normally considered. 

The results derived here concern properties of 8 th order symmetric tensors, and the number of 
constants they possess under transverse isotropy and isotropy, subject to the indicial symmetries 

We begin in Section [2] with a summary of the main results. Although the focus of the paper 
deals with transformations caused by rotation, it will prove useful to first consider invariance under 
reflection. Section[3]introduces the notion of invariance under reflection about a plane, and quadratic 
invariants of the elastic moduli are derived in Section H] for transformations under one and two 
reflections. Transformations under rotation about an axis is then considered in Section [5] where the 
results are proved. Finally, issues of consistency and completeness are discussed in Section El 

A note on notation: The vector triad e^, i — 1, 2, 3 is an ortho normal basis in 3-dimensions. The 
summation convention on repeated suffices is assumed. 



2 Summary of the principal results 

The results of the paper are summarized in the form of two theorems, with the proofs given in the 
subsequent Sections. 

Define the 21-vector of elastic moduli 

C = (en C22 C33 C23 C13 C12 C44 C55 C66 C14 C25 C36 C34 C15 C26 C24 C35 C\q C$q C46 C45) , (12) 

then we have the following: 

Theorem 1 All quadratic invariants of c under SO(3) are of the form c l Fc where the 21x21 
symmetric matrix F is of the form 

jp _ ( ^9x9 0g x i2 
\0l2x9 -&L2xl2 

with symmetric matrices A and B defined by 12 and 7 independent elements, respectively. These 
matrices have the form 



I hi 


/l2 


fl2 


fu 


fl5 


fl5 


fll 


/l8 


fl8 


\ 




hi 


fl2 


fl5 


fu 


fib 


fl8 


/l7 


fl8 








fll 


fl5 


fl5 


fu 


fl8 


/l8 


hi 










fte 




fi5 


fi7 


^48 


fl8 












^44 


hb 


^48 


At 


f±8 














f\l 


^48 


^48 


fn 






s 


Y 


M 








^78 


fl8 


















fn 


fl8 




V 
















fn 


) 



(14) 
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B 
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V 



and the 19 distinct elements 









bu 








bu 








&110 








b n 








bu 








bu 








biw 







b n 








bu 








bu 








h 10 






bu 








hi 








&410 














644 








647 








&410 













644 








647 








b± 10 












644 




644 





644 


&410 







&410 





b± 10 


S 


Y 


M 












frioio 




frioio 





frioio 



(15) 



fill fl2i fui /l5) /l7) /l8, fiii ^45; ^47; ^48) fill fl%i ^11) ^14, &1 ^44, &47) ^4,10) ^10, 10) 

are arbitrary as long as they satisfy the 12 linearly independent conditions 



2/11 


+ 2/ 12 -4/ 15 + / 44 


— /77 


= 0, 


(16a) 


2/11 


+ 2/ 12 — 2/ 18 — / 44 


-hi 


= 0, 


(16b) 




2/12 — flA 


-hi 


= 0, 


(16c) 




fl4 + fl5 — ^45 


— hs 


= 0, 


(16d) 




/l7 + A 8 — ^48 


— hs 


= 0, 


(16e) 




2/44 — 2/45 


-bn 


= 0, 


(16f) 




2/l4 — 2/i5 


+ &14 


= 0, 


(16g) 




2/47 — 2/48 " 


- h 10 


= 0, 


(16h) 




4/n-2/ 16 -2/ 18 


- ^44 


= 0, 


(16i) 




4/i2- 2/15 -2/ 18 


+ &4 7 


= 0, 


(16j) 




2/l7 — 2/i 8 " 


f &410 


= 0, 


(16k) 




2/77 — 2/78 - 


&1010 


= 0. 


(161) 



The number of linearly independent quadratic invariants under SO (3) is seven. 

An immediate corollary of Theorem [J is that the seven invariants derived by lAhmadl ( 20021 ). 
based partly on the work of iTingl f 19871 ). are linearly independent and complete. The second result 
is 

Theorem 2 All quadratic invariants of c under SO (2) are of the form c l Fc where the 21x21 
symmetric matrix F is 

F = F (1) +F (2) ; 

and F^ define subspaces of dimension 29 and 6, respectively, and 



'A 



9x9 







9x12 







_p(2) 



,0 



12x9 



B 



(1) 
12x12/ 



9x9 



D* B^ 

■-^12x9 - D 12xl2/ 




The symmetric matrices A, B^ and B^ and the matrix D are defined by 17, 27, 3 and 3 inde- 
pendent elements, respectively. 
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With the axis of rotation e 3 , 



A 



/ f 
' Jll 


f 

Jl2 


Jl3 


J14 


J15 


f 

Jl6 


f 

J17 


Jl8 


J19 


\ 




f 

III 


f 

Jl3 


/ 

J15 


/ 

J14 


f 

Jl6 


/ 

Jl8 


f 

In 


f 

J19 








f 

J33 


f 

J34 


f 

J34 


f 

J36 


f 

J37 


J37 


f 

J39 










f 

J44 


f 

J45 


f 

J46 


J47 


J48 


J49 












J 44 


J 46 


J 48 


J 47 


J 49 














^66 


^67 


fe7 


/69 






s 


1" 


M 






^77 




/79 


















fn 


/79 




V 
















/99 


/ 



(19) 



11 









b u 








fel7 








h 10 
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fell 







617 






bu 








61 10 
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&36 








he 








^312 








fe 44 






647 








hlO 















&55 
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hn 
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bm 








bei2 
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hll 






















fe 44 







bi 10 





















bm 







bei2 






Y 
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&1010 




























frioio 



























b\2V2 


) 



(20) 



V 

The elements fs3, ^34, ^37 0/ A and 6 44 of B^ are arbitrary, while the remaining 40 distinct 
elements in A and are arbitrary as long as they satisfy the 15 linearly independent conditions 



2/11 


+ 2/ 12 -4/ 16 + / 66 


— f99 


= 0, 


(21a 


2/11 


+ 2/12 — 2/ig — fee 


— fe9 


= 0, 


(21b 




2/13 — fze 


— f'39 


= 0, 


(21c 




/14 + fi5 — fie 


~ fi9 


= 0, 


(21d 




fx 7 + h s — fe 7 


— /79 


= 0, 


(21e 




2/44 — 2/45 


- &33 


= 0, 


(21f 




2/i4 — 2/15 


+ ^36 


= 0, 


(21g 




2/47 — 2/48 - 


~ h 12 


= 0, 


(21h 




4/ii-2/i 6 -2/i9 


- ^66 


= 0, 


(21i 




4/ 12 -2/i 6 -2/ 19 


+ ^69 


= 0, 


(21j 




2/17 — 2/i 8 - 


t" &612 


= 0, 


(21k 




2/77 ~ 2/ 78 — 


&1212 


= 0, 


(211 




6l 4 — &4 7 " 


I" ^410 


= 0, 


(21m 




fell - 26i 7 + 655 - 


frio 10 


= 0, 


(21n 




6ll + &110 - ^55 " 


1- hu 


= 0. 


(21o 



NORRIS 



Quadratic invariants 



6 



Also, 
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-di3 
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-di 12 
















































\ 



D 



J 

where the 6 elements 615, ^45, d\3, di t \2, ^4, 12 ore arbitrary. 

The number of linearly independent quadratic invariants under SO (2) is thirty five. 



(22) 



(23) 



3 Invariants under reflection about a plane 

3.1 Linear and quadratic invariants under reflection 

In order to fix ideas consider a tensor of order one, a vector v = v^i. Invariants are defined 
in relation to a transformation, usually associated with some material symmetry. The simplest 
transformation is that of reflection about a plane orthogonal to a direction e, which we denot e 
R(e). This transformation is associated with monoclinic symmetry ( Cowin and Mehrabadil . fl995h . 
The action of R(e) on v is defined by 

R(e)y = v - 2(v ■ e)e. (24) 

Thus all v orthogonal to e are unchanged, or invariant under R(e). The linear invariant of an 
arbitrary vector may be defined as its projection onto the invariant subspace, i.e. v — (v • e)e. The 
number of linear invariants is the dimension of the invariant subspace, in this case 2. 

Quadratic invariants are defined by quadratic forms, which for a vector require a symmetric 
second order tensor, say F = F . Let 

0(v) = v-F-v, (25) 
then we seek F which leave <j) unchanged under the action of R(e) on v: 

= R(e)4> where R(e)<f> = <f>(R(e)v). (26) 
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Since 

R(e)(p = [v - 2(v ■ e)e] • F • [v — 2(v ■ e)e] = - 4(v • e) [v - (v ■ e)e] • F • e, (27) 

we see that R(e)<fi = <p under two circumstances: (i) for all F such that F • e = 0, which defines a 
three dimensional subspace of second order symmetric tensors; (ii) for all F of the form F = fie ® e, 
a one dimensional subspace. It may be checked that all quadratic forms that leave <ft fixed must be 
combinations of these, and hence the number of quadratic invariants for vectors is 3+1=4. 

As we deal with tensors of higher order it is simpler to work with the components relative to 
the basis e i; i — 1,2, 3. To be specific, we consider e = e 3 , then the 2 linear invariants are V\ and 
V2- Similarly, the 4 quadratic invariants are defined by <f> of the form (125]) where /13 = j'23 = but 
F is otherwise arbitrary. In summary, 

fvi\ ffn fu 0\ 

v = Uj , F = /ia /aa , n = 1 under R(e 3 ), (28) 
W \0 f 33 ) 

where n — 1 indicates that these are the linear and quadratic invariants for tensors of order n. 

Note that we do not present the quadratic invariants as, for instance, the set {v\, v\, v 1V2, v 3 }, 
but use notions from linear algebra which are natural for quadratic forms. 



3.1.1 n=2 



A second order symmetric tensor defines a unique quadratic form, through eq. 
analysis for n = 1 we can immediately identify the linear invariants under R(e 3 
symmetric tensor V 



e,- as v 1 i,v 2 2,v 12 ,v 33 . 



Quadratic forms on V are 



where the elements of F satisfy 



VijFijkiVki, 



231) . Based on the 
any second order 



(29) 



ijkl 



jikli 



ijkl 



F, 



klij ■ 



(30) 



These are the same relations that define Ela, see eq. (Q. Thus, the symmetries of F^i under R(e 3 ) 
follow by standard arguments: all el ements with index 3 occurring an odd number of times must 
vanish ( Cowin and Mehrabadil . fl995l ) . In summary, 



V 



v n v 12 
v±2 v 22 

t> 33 , 



ffn 


/12 


/l3 










fu 


/22 


/23 








/i6 


fl3 


/23 


/33 








/36 
















/46 














/55 





U16 


/26 


/36 











n = 2 under R(e 3 ). (31) 



The thirteen elements of F define the subspace of Ela invariant under the single reflection 
R(e 3 ). That is, there are 13 linear invariants for fourth order symmetric elasticity tensors. Since 
these tensors are defined by the quadratic form (energy) acting on second order symmetric tensors 
(strain) the 13 linear invariants of n = 4 correspond to the quadratic invariants of n = 2, in the 
same way that the linear invariants for n = 2 correspond to the quadratic invariants of n = 1, eqs. 
and (l3~T]) i . It should be clear that the same equivalence holds for arbitrary tensor order n, 



and the result may be stated as follows: 



Lemma 1 The quadratic invariants of a tensor of order n are defined by the linear invariants of 
the corresponding symmetric tensor of order 2n. 
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3.2 Remarks on reflections and rotations 

We are concerned with invariants under the group of proper orthogonal transformations, SO (3), and 
the group of rotations about an axis, SO (2). The orthogonality requires that right-handed triads of 
coordinate axes remain right-handed under transformation, or equivalently that the determinants of 
the 3x3 transformation matrices are +1. But three dimensional inversion is also acceptable by virtue 
of the fact that the tensor we are concerned with is of even order, viz. 8. Therefore, 5*0(3) should 
be considered extended, to include the group J 3 of transformations 1,2,3 — > —1,2,3, etc. Thus 
0(3) = 50(3) U 1%. Similarly, 0(2) = 50(2) U I2 where I2 is the group of transformations formed 
from 1,2 — > —1,2; 1,2 — > 1,-2. Another way to understand why we can replace 50(2) — » 0(2) 
is to consider reflection about a plane, say -R(e 3 ), followed by inversion. This is equal to rotation 
about the normal to the plane by ir, and hence 0(2) is equivalent to SO(2)Ui?(e 3 ). We will find 
this particularly useful as a starting point later, and it motivates consideration of the invariants 
under the action of i?(e 3 ) first. 

We are now ready to determine the quadratic invariants of Ela, starting with invariants under 
reflection. 

4 Invariants of elastic moduli under reflection 

We examine Ela under reflection first about a single plane, and then about two. Based on the Voigt 
notation for indexing elastic moduli with 21 elements, the indices for F of ( fTOj) and ( fTTI) as a square 
symmetric matrix run from 1, 2, . . . , 21, according to the following 

i = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 
IJ = 11 22 33 23 13 12 44 55 66 14 25 36 34 15 26 24 35 16 56 46 45 
where IJ are the Voigt indices, which are capital suffices taking the values 1, 2, . . . , 6 by 

1 = 1,2,3,4,5,6 ij = 11,22,33,23,31,12. (32) 

For example, /1317 = /343s = /33233313 using the two, four and eight index notations. In this way 
the elements fijkipqrs are represented uniquely by the ^(21 x 22) = 231 elements /y. 

Consider first the reflection R(e 3 ) about the X1X2 plane as a symmetry plane. Since we are 
dealing with a tensor of even order (8), it follows that those elements vanish that have indices 
ijklpqrs where 3 occurs an odd number of times. There are 13 x 8 = 104 elements in this category: 

f j3k3 = 0, j 3 G {1, . . . , 9, 12, 15, 18, 21}, h e {10, 11, 13, 14, 16, 17, 19, 20}. (33) 

That leaves 231 — 104 = 127 non-zero elements that remain under reflection about a symmetry 
plane, or monoclinic symmetry. 

Next consider reflection about two planes with perpendicular normals, say -R(e 3 ) and R{g\). As 
before, elements that have indices ijklpqrs where 1 occurs an odd number of times must vanish. 
The 104 elements with 1 occurring an odd number of times are 

f jlkl =0, jt E {1,..., 9, 10, 13, 16, 19}, fcie {11, 12, 14, 15, 17, 18, 20, 21}. (34) 

But some of these coincide with members of (I3"3"|) : the elements indicated by the indices in boldface, 
plus the elements defined by the remaining indices. The elements that are distinct from ( 133]) are 
the 9 x 4 = 36: 

/JA = 0, jiG{l,...,9}, ^€{12,15,18,21}, (35) 

and the 4 x 4 = 16: 

f hh = 0, j x G {10, 13, 16, 19}, k e {11, 14, 17, 20}, 
In summary, F is of the form 

F _( ^9x9 0g x i2 
\0l2x9 -&L2xl2 



(36) 
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where A and B are symmetric. The matrix A is full, indicating 9 x 10/2 
and letting the indices of B run from 1 to 12, we have 



45 independent elements, 



0, (ij) G {{(3, 6, 9, 12) x (1, 2, 4, 5, 7, 8, 10, 11)} © {(1, 4, 7, 10) x (2, 5, 8, 11)}}. (38) 

4 x 12 = 30 for a total of 75 independent elements in F. In fact, B 



Therefore, B has 12 x 13/2 
has banded structure 
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^12 12 
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(39) 



run from 1, 2, . . . , 12, according to the 



Note the indices for B as a square symmetric matrix bi d 
following 

B i = 123456789 10 11 12 
Voigt IJ = 14 25 36 34 15 26 24 35 16 56 46 45 
F(21x21) i = 10 11 12 13 14 15 16 17 18 19 20 21 

The 8 th order tensor F of eqs. ( !37l) and ( l39l) is unchanged under reflection about a third plane 
orthogonal to the others, i.e. R(e 2 ). In order to see this, note that R(e 2 ) implies that all elements 
with index 2 occurring an odd number of times must vanish. However, since the tensor is of even 
order, and we have eliminated elements with 1 and 3 occurring odd numbers of times, the tensor 
is automatically invariant under R(e 2 ). This effect is well kno wn in elasticity: that a materi al with 
two orthogonal planes of symmetry automatically has a third f lCowin and Mehrabadil . Il995l ). 

In summary, 

Lemma 2 All quadratic invariants of c under reflection about a plane are of the form c l Fc where 
the 21 x 21 symmetric matrix F is 



-9x9 

t 

12x9 



-Dg x i2 
B\2xl2 



(40) 



A = A 1 is full and B = B l and D are of banded form, 



( b u 



B 



bi2 

b 2 2 






bu 


bi5 





bn 


bis 





bi 10 


bin 








b 2 4 


&25 





b 2 7 


&28 





b 2 10 


&211 





b 3 3 








b 3 6 








&39 








&312 




b M 


b 4 5 





647 


b 4 s 





^410 


b 4 n 









b 5 5 





b 5 7 


^58 





bsio 


hn 











bee 








bm 








be 12 










b 77 


b 7 s 





^710 


b 7 n 















^88 





bsio 


&8ii 

















bgg 








bgu 
















^1010 


bion 





















fen 11 























b\2!2 



(41) 
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and 



D 



1 





di3 
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dig 








"1 12 


\ 








<^23 








"26 








G?29 








"2 12 










<^33 








"36 








C?3g 








G?3 12 










d^3 








C?46 








C?49 








dti 12 










d 5 3 








^56 








d§9 








<^512 










d 63 








^66 








d&9 








d§Y2 










d73 








d76 








d79 








0^712 










dg3 








dge 








ds9 








0^812 




V o 





d93 








"96 








dgg 








"912 


/ 



(42) 



where the second index of D is the same as that for B, i.e. runs from 1,2, ... , 12. The number of 
linearly independent quadratic invariants is 127. 

Under reflection about two orthogonal planes the number of linearly independent quadratic in- 
variants reduces to 75, and F then has the form with D = and the 16 elements in bold in 
(f^iy are zero, i.e. B reduces to B^ of eq. ( dZ| ). 



Thus, invariance under reflection about one and then two orthogonal planes reduces the number 
of independent elements in F from 231 to 127 to 75. The case of a single plane is monoclinic 
symmetry, and two (or three) corresponds to orthorhombic symmetry. 

Reflection symmetry is an i mportant tool in developing the structure of anisotropic elastic tensors 
( Cowin and Mehrabadil . Il995l ). In fact, all eig ht fundamental symme tries can be cast in terms of 



repeated application of the reflection operator (jChadwick et al.l . 200 ll ). The definition of transverse 
isotropy implies that every plane containing the axis of symmetry is a symmetry plane, in addition 
to the plane orthogonal to the axis. However, this equivalence does not hold in the present problem, 
dealing with the form of the 8 th order tensor F invariant under SO (2). Consider rotation about e 3 
by 7r, which transforms {ei,e 2 } — > {e 2 , — ei}. Since elements with odd numbered occurrences of 3 
have already been eliminated by the requirement of invariance under -R(e 3 ), it therefore eliminates 
terms with 1 (or 2) occurring an odd number of times. Thus, C3312 = for TI elasticity, but it does 
not eliminate elements such as /33122233 from F. Hence SO (2) invariance for F is not the same as 
invariance under R(e) for all e _L 63. 

In other words, SO (2) invariance of F is not a subspace of orthorhombic symmetry, but it is a 
subspace of monoclinic symmetry. This allows us to start the search for SO (2) invariance with only 
127 of the 231 elements. Furthermore, SO (3) invariance of F is a proper subspace of orthorhombic 
symmetry, meaning it can be sought starting from the 75 element form of F in ( 1651 . 

First we need to introduce the rotation matrix for Ela. 



5 Invariants under rotation about an axis 

We commence by casting the problem in general terms, with some simpler examples before consid- 
ering the problem for the 8 th order tensor. 



5.1 General theory 

Let c be a vector, not restricted to the elastic modulus vector, and consider the quadratic form 

= c'Fc, (43) 

where F is symmetric^] 

(44) 



F\ 



2 If F = M + S where M = M t and S = -S l then defined in (gHJ is 
skew-symmetric part of F is irrelevant and can be set to zero. 



c*Fc independent of S. The 
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Under rotation through angle 9 about axis p, a 3-vector, 

c^d = Qc, (45) 
where Q = Q(9,p) is the rotation. Accordingly, 

_> fttf) = c'Q'FQc. (46) 
We want 4>'(9) to be independent of 9, that is, its derivative should vanish for all 9. The derivative 



is 



0' = c'Q'FQc + c'Q'FQc, (47) 

where the dot indicates differentiatio n with respect to 9. By an application of standard Lie group 
and Lie algebra theory (jFeganl . Il99ll ). Q{9) is a one parameter subgroup of SO(N) where N is the 
vector length. Q is generated by P = — P 4 , an element from its Lie algebra of skew-symmetric 
matrices. Thus, Q{9) = exp9P, from which the derivative is 

Q = PQ. (48) 



The particular form of P for elasticity tensors (N = 21) was derived by iNorril ( 2007 ). 
Equations (}4"Tj) and f|48|) imply 

4/ = c l (Q'FPQ - Q'PFQ) c. (49) 

Note that the matrix (Q l FPQ — Q t PFQ) is symmetric, hence <f)'{9) = for all possible c implies 
the SO (2) invariance condition: 

Q\FP - PF)Q = for all 9. (50) 

Hence we have 

Lemma 3 The quadratic form <fi = c t Fc is an SO (2) invariant if and only if F commutes with the 
generator P: 

FP-PF = 0. (51) 

Equation (|5T|) is a necessary condition since (|50|) must hold at 9 = where Q(0) = 1. The sufficiency 
of flBT]) is clear, hence eq. (15T|) is the condition that F must satisfy. 

The question of quadratic invariants under SO (2) therefore reduces to finding all symmetric 
F which commute with the skew rotation generator P, or equivalently, sym(FP) vanishes. We 
next demonstrate this approach by application to increasingly higher order tensors, associated with 
quadratic forms on tensors of order n = 1, 2 and ultimately the desired n = 4. 

5.2 Quadratic forms on vectors: n=l 

In this case the generator is the standard and well known skew symmetric second order tensor 
defined by the axis of rotation p, 

Pij = -tijkPk- (52) 

Let p = e 3 so that 




and let 

fll fl2 fl3 

f\2 /22 /23 I • (54) 



,/l3 ^23 /: 



33 , 
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Then 

( — 2/12 fu — fzi —fi?\ 

fu - / 22 2/ 12 / 13 . (55) 

-/ 2 3 /13 / 

Setting each of the elements to zero gives four conditions: fu — fii = /12 = /13 = fzz = 0, so that 
the most general F is of the form 

'a N 

a I , (56) 
,0 b, 

for a, b 7^ 0. 

The associated invariants are 

a = c\ + c 2 , b = C3. (57) 

For arbitrary direction p these are 

a = C * C -(pV) 2 , fo = (^c) 2 , (58) 

corresponding to 

F a = -P\ F b = I + P 2 . (59) 

Any F comprised of even powers of P will obviously commute with P. For n = 1 there are only 
two independent F of this form, i.e. F — I, P 2 , since P satisfies the characteristic equation 

p 3 + p = 0, (60) 

and hence P 2+2m = (-l) m P 2 . 



5.3 Quadratic forms on tensors: n=2 

The relation between the second order symmetric tensor C and the 6-vector c is 

> c = (C u C 22 C 33 V2C 23 V2C 13 V2C 12 y 
P is now a 6x6 skew symmetric matrix, 

y/2(Y - Z l 

y/2{Z — Y l ) X - X 1 





(Cu 


Cl2 


c 




C12 


c 22 


c 




VC13 


C*23 


C 



(61) 



P 



(62) 



Y 



(63) 



where 

'0 p 2 0\ /0 
p 3 , Z=\Q 

ft O/ \P3 

The m atrix P was first derived by iMehrabadi et al.l ( 19951 ). and the present form is due to iNorris 
( 120071 ) . We will find this useful when we consider the anal ogous P for n = 4 later. We note that the 





characteristic equation of P is now of fifth (2n+l) order (IMehrabadi et all Il995l : iNorrisl . 120071 ) 



P(P 2 + I){P 2 + 41) = 0. 



(64) 
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Let p = e 3 then 



( 














-V2\ 

















V2 








(J 























1 














-1 








\V2 


-V2 











o J 



(65) 



and the SO (2) commutator is 

/— 2/i6 fw — I26 

2/26 

FP-PF 



V2 







775/15 ~ /46 — 77^/14 — /: 



75/25 + fi6 



y^J14 ~ ,7 56 
— 75/24 + /56 

~72^ 34 
75 (/55 — /44) 



5' 



1' 



M 



-V2f, 



15 



/ll — /l2 ~ /66 \ 
/l2 — /22 + /66 
/l3 — /23 
/l4 — /24 + 75/56 
/l5 — /25 — :75/46 



V 



V2> 
2/l6 — 2/26 



(66) 



Setting all elements to zero and solving the resulting linear equations, we find that F has five 
independent elements (and the expected form of a TI tensor of elastic moduli), 



/ In fl2 I13 

fu /11 /13 

/l3 /l3 /33 





\ 

Thus, there are five quadratic invariants, 







/44 











/44 













(67) 



/11 ^ /12/ 



>n — C*ii + C22 + 2Cf 2 , 0i2 — 

013 = (Cll + ^22)^33, 044 



CllC22 



c 2 

r<2 1 ^2 

°13 °23) 



12; 



^33 



r 2 



(68) 



corresponding to f u , /12, /33, /13, and / 44 , respectively. 
The linear invariants of C under SO (2) are 



Ax = C 



11 



22- 



Ao — C; 



33- 



(69) 



The five quadratic invariants can be recast as three defined by the linear invariants and two new 
ones, e.g., 



>1 = \j, 02 



X 2 
A 2; 



A1A2, 



C11C22 — c 



12; 



Cl3 + ^23- 



(70) 



In keeping with the statement of Lemma [U we note that the five independent elements of F 
correspond to the linear invariants of Ela under SO (2). We next consider the quadratic invariants 
of Ela under SO (2). 
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(71) 
(72) 



5.4 Quadratic forms on the elastic moduli: n=4 

The fundamental quadratic form is now expressed in alternative ways, 

(f) = jFc = c l Fc, 

where c is a 21-vector related to the 21-vector of moduli c in (|12[) by 

c = Tc, 
and 

T = diag (1 1 1 V2 V2 \/2 2 2 2 2 2 2 2 2 2 2 2 2 V8 \/8 \/8). (73) 

We introduce F in (ITT]) to remove the \/2 factors from the final expressions. The two are simply 
related by 

F = TFT. (74) 
Results below are given for the elements of F which do not have the \/2 factors. 

5.4.1 The rotation matrix 

The 21x21 skew symmetric generator is ( Norrisl . [20071 ) 

P = R-R\ 



(75) 



where 



R 



(o 











2Y 

















-V2Y 

























2N 


-V2Y 





-V2Z 








X 





-V2X 





V2N 


2N 








X 





2Z 








X 








V2X 







-V2Z 


-V2X 


V2X 





-X 



(76) 



with X, Y and Z as before in eq. (1B"3"|) . and 



N 



( V\ o 
p 2 
p 3 , 



Note that the characteristic equations for P is 

P(P 2 + I)(P 2 + 4/)(P 2 + 9I)(P 2 + 16/) = 
The conditions (|5"T|) are that the commutator of F and P vanish, 

FP-PF = 0. 



0. 



(77) 



(78 



(79) 



This may be rewritten in a form involving the elements of F of (174]) . with the intent of removing 
the factors of \/2. The resulting equation is 



PF + FP l 



0. 



where P = TPT' 1 . 



(80) 



The modified generator P is no longer skew symmetric, but it still satisfies the characteristic equa- 
tion (ITS"]) , and, like P, can be used to compute the rotation Q = exp(6P), or its modified form, 
Q = exp(6>P) = TQT" 1 , as a polynomial of degree 8 in P (jNorrisl . 120071 ). With p = e 3 , P has 18 
non-zero elements, hence P has 36 non-zero elements. 
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5.4.2 Quarter turn conditions on the elements of F 

Before applying the SO (2) condition (]80[) to F, we perform a preliminary simplification by invoking 
invariance under a quarter turn about the axis of rotation. The transformed elements Fijki pqrs should 
be unaltered under the interchanges of indices 1,2 — > 2,-1. This implies a total of 52 relations 
among the 127 elements of F, which we split into two categories depending on whether the index 1 
occurs an even or an odd number of times. In the former category are the following 31 identities, 



fll 


— ^22, 


/l3 


= /2 3, 


/l4 


— /25, 


7l5 


— /2 4, 


/l6 


= /26> 


fl7 


= /2 8, 


/l8 


— /2 7, 


fl9 


= /2 9, 


hi 


= /35j 


fz7 


= /38, 


/44 


= /55, 


/4 6 


= /56, 


fi7 


= /58j 


/4 8 


— /57, 


/4 9 


= /59, 


/6 7 


= /68, 


J77 


= ^88, 


A 9 


= ^89, 










fell 


= b 22 , 


fel4 


= fe28, 


fel7 


= fe2 5, 


fel 10 


= 6211 


he 


= ^39, 


fe 44 


= fe88, 


fe 47 


= fe58, 


fe410 


= fesn 




= b 7 7, 


6511 


= feno, 


fee 6 


= fe99, 


fe612 


= fe912 


1010 


= fell 11, 















(81) 



and in the second group we have the following 21 connections, 



fel 5 


+ fe 2 7 = 


0. 


fel 8 


+ fe 2 4 = 


0. 


fel 11 


+ fe2 10 = 


: 0, 


&4 5 


+ b 78 = 


0. 


64 11 


+ fes 10 = 


= 0, 


fesio 


+ &711 = 




di 3 


+ d 23 = 


0. 


d*3 


+ d 53 = 


0. 


4 3 


+ 43 = 


0. 


di e 


+ d 29 = 


0. 


dig 


+ d 26 = 


0. 


^3 6 


+ 49 = 


0. 


C?46 


+ d 59 = 


0. 


C?4 9 


+ d 56 = 


0. 


dee 


+ d^g = 


0. 


<i 76 


+ d 89 = 


0. 


0?79 


+ d S 6 = 


0. 


dg§ 


+ dgg = 


0. 


^112 


+ d 2 i2 = 


= 0, 


^412 


+ 4 12 = 


= 0, 


d-j\ 2 


+ d 8 12 = 


= 0. 



(82) 



Therefore, A, B and D have 45-18=27, 46-19=27 and 36-15=21 independent elements, respec- 
tively, for a total of 127- 52= 75 independent elements in F . A has the form given in (1191) . while B 
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and D are of the form 
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(83) 



5.4.3 SO(2) conditions on the elements of F 

We are now ready to apply the SO (2) invariance condition flHUj) to the simplified form of F with 75 
independent elements, where the matrix P defined by (1761) with p = e3 has 36 nonzero elements. 
The system of equations generated is large but relatively straightforward to solve, particularly with 
the aid of an electronic computer. We simply state the results, which the reader may check. 

We find that the matrix in flHUj) vanishes if and only if 40 additional relations among the 75 
independent elements are met. These may be split into two sets, of 15 and 25 respectively. The first 
set of 15 relations (T2T1) involve only those elements of F in which the index 1 (and hence 2) occurs 
an even number of times, and therefore they act only on the elements of A 9x9 and the non-zero 

elements of B^xu °f (120]). The second set of 25 relations involve only elements of F with the index 
1 (and 2) occurring an odd number of times: 

&i,2 = 0, 64,8 = 0, 65,7 = 0, 610,11 = 0, 

&1.5 - &1,H = 6l,5 + & 5,10 = 0, 6i >8 + &4,5 - &4.11 = 0, 



4,6 


= 0, 


4,9 


= 0, 


4,3 


= 0, 


4,6 = 


0, 


4,12 


= 0, 


4,3 


= 0, 


4,3 


= 0, 


4,6 — 


0, 


4,12 


= 0, 


4,12 


= 0, 


4,3 


= 0, 


4,6 — 


0, 


4,12 


= 0, 


4,3 — 4,6 


= 0, 


4,3 ~~ 4,9 


= 0, 


4,12 — 4,6 — 


0, 


— 4,9 


= 0, 


4,12 + 4,3 


= 0. 











(85) 
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Thus, 
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(86) 

and D is given by (|23|) . The terms in bold in (1861) are essentially decoupled from the others, and we 
therefore split F according to eq. (ITT)) in order to emphasize the disjoint nature of the subspaces 
generated by the elements with index 1 occurring an even and an odd number of times. 

The partition (ITTj) also allows us to easily determine the dimensionality of F. Thus, F^ clearly 
has 6 independent elements. There are 44 distinct elements in F^\ and all but 4 of these, ^33, ^34, 
^37 and 644, occur in the relations 02 These 15 conditions are of rank 15, indicating that the 
conditions are linearly independent. This can be verified by noting that each of the 15 equations 
involves at least one element (the final one in the left member) not contained in the other 14 
equations. Hence, there are a total of 15 constraints on the 50 elements of F, which therefore has 
50-15 = 35 independent elements. This is the dimension of F, and the number of independent 
quadratic forms on Ela. 



5.5 SO(3) conditions on the elements of F 

There are various ways to deduce the SO (3) form of F using the results for SO (2) (actually 0(2)). As 
discussed before, the required group of transformations is simply 0(3). We also note the following 
as a consequence of Theorem [2] and Lemma [2j 

Lemma 4 All quadratic invariants of c under 0(2)UR(ej_) where e± is perpendicular to the axis 
of rotation are of the form c l F^c where the 21x21 symmetric matrix F^ is given by eq. {\18§ i . 
The number of linearly independent quadratic invariants under SO(2)UR(ej_) is 29. 

We can therefore start with F^- 1 ' of Theorem [2], which corresponds to 0(2)Ui2(ei). The next 
step is to consider quarter turns about ei and e2, which reduces A and B to the forms (114)1 and (115)) . 
respectively. The 12 distinct elements in A and 7 in B are related by the 15 conditions of Theorem 
|2]in eqs. ( l2~Tj) . We find that only 12 of these are linearly independent, or in other words, the system 
of equations is rank(12). A linearly independent system can be obtained by, for instance, ignoring 
the final three conditions in (1211) . to give a system of 12 equations on the 19 elements, i.e. eqs. ( fl6l) . 

In summary, there are 12 relations between the 19 elements of F. Hence there are 7 quadratic 
invariants under SO (3). 

It remains to discuss these results with respect to the invariants proposed by iTingi (jl987l ) and 



Ahmadl ((20021) 



6 Consistency and completeness 

A given quadratic form can be checked to see if it is consistent with one of the invariant forms 
defined in Theorems [D and El We will describe how to do this and discuss the consistency of the 
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quadratic invariants proposed by Ting and Ahmad. We note in passing that A±, A 2 and Bi, . . . , B4 
of eqs. (ED) and (|2j) are obviously consistent with SO (3) invariance, but some of the SO (2) invariants 
proposed by Ting and Ahmad are not immediately obvious. We also discuss the completeness of 
Ting and Ahmad's quadratic invariants. 



6.1 SO(3) 

6.1.1 Consistency 

Given a quadratic form in c, for instance A\ of eq. (pQ), define a 21x21 symmetric matrix F 
according to 

f - 1 ^ (Z7) 

fij ~ 2dc~b\/ (87) 

where q is the i— th component of the modulus vector c of ( 1T21) . The quadratic form is an invariant 
under SO(3) iff (i) the elements of F have the form as defined by eqs. 013p - (|15p . and (ii) they satisfy 
the 12 equations (fTBI) . 

It may be confirmed that Ahmad's seven quadratic forms {A\, A\, A\A 2 , B±, B 2 , B 3 , 5 4 } are 
SO (3) invariants. 



6.1.2 Completeness 

Assuming a quadratic form is consistent with SO (3) invariance, define the vector of the 19 distinct 
elements in F of (fl4l) and ( fl5j) . specifically 

U = «19xl = (fll fl2 /l4 fl5 fll fl8 fiA fib fil fi8 fll fl8 ^11 ^14 ^1 10 ^44 ^47 ^4 10 ^10 lo) • (88) 

Now suppose we have a set of different quadratic forms, fa, . . . , 0jv, where N < 7 (we do not need 
to consider N > 7 since there can be no more than 7 linearly independent forms). Let Ui be the 
19- vector for fa, i — 1, . . . , N, and define the matrix 

M = M 19xN = ( Ul u 2 ...u N ). (89) 

Then the quadratic forms are linearly independent iff rank(M) = N. By definition, a set of 7 
quadratic forms is complete if they are linearly independent. 

It may be checked that Ahmad's seven quadratic forms define a matrix M of rank 7, and are 
therefore complete. 



6.2 SO(2) 

Checking a given quadratic form for consistency is analogous to the procedure described for SO (3). 
Thus, first compute F according to flH7|) . then check that (i) the elements of F have the form as 
defined by eqs. f|T3|) - f|T5|) . and (ii) they satisfy the 15 equations ff2~Tj) . 

It may be confirmed that Ahmad's 17 quadratic forms Ei, . . . , E 17 in eq. ([7]) and the new 
quadratic forms E i8 , E 19 , E 20 are all SO (2) invariants, as are the 15 quadratic forms defined by the 
five linear invariants L\, . . . , L5 of eq. ([6]). 

These 35 quadratic forms are also complete. In order to see this define v = v 50xl as the vector 
of the 50 distinct elements in F of ffT9l) - (123]) . specifically 

v — {fll fl2 fl3 flA fib fw fll fl8 fl9 f33 fzi ^36 f'il ^39 ^44 fi5 fi6 

fvi fm ^49 fm f%i fm fii fis fid fm &n ^14 ^15 &17 &18 &110 ^33 

&36 ^312 &44 &45 ^47 ^4 10 ^55 ^511 ^66 ^69 ^612 &10 10 ^12 12 ^13 ^112 ^412) • (90) 
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Let M be the matrix formed from the w-vectors of the 35 quadratic forms. Then it may be checked 
that rank(M) = 35, indicating that they complete set. 

Finally, we note that 

-jF^c = -d 412 E 15 - d 13 E 16 + d 112 E 17 + b 15 E 18 + (6 18 + 645)^19 - b 18 E 20 . (91) 

This illustrates how the partition of F in eq. (ITT]) splits the set of invariants into distinct subsets, 
one associated with elements of F that have index 1 (and 2) an even number of times, which is a 
29 dimensional subspace. The other is the six dimensional subspace associated with the fact that 
SO (2) is not a subspace of orthorhombic symmetry for this type of 8 th order tensor. 



Appendix 

A Proof of equation (3D 

Define the totally symmetric and the asymmetric parts of the elasticity tensor as 



Cijli — -^(Cij k i + C ik ji + Cujk), Cijli ~ g(2Cyfc* — Ci^i — Cujk), (A.l, 



respectively. These partition the elastic moduli 

'ijkl "T" ^ijkV 



Cij ki — C\ m + C- • ^ , ( A . 2 ) 



and are orthogonal in the sense that 

.n. ... — r 

ijkl ijkl ~ ijkl ijkl' 



n n — _i_ r*( a ) l a r>\ 

^ijkl^ijkl — ^ijki^ijkl u ijkl u iikl- K^-d) 



Hence, 



and therefore 



2CijkiCikji — CijkiiCikji + Cujk) 

= Cijki(3C^l t — Cijki) 

— r* / , 9r | ( s ) n( a ) \ 

— ^ijkii^ijki - ^ijki) 

_ o/i( s ) r*^ _ r*(°) r*^ 

— ^ijkV^ijkl ^ijkl^ijkl 

= ICijuCijki - %Cfj U C^ u , (A.4) 
B 5 = B 1 - 3 -C^ kl C% (A.5) 



While it is easy to verify that 



Cl-l = \Dij, (A.6) 



where 



F>ij — Cijkk ~ Cikjk-, (A.7) 
it is not as obvious that C^ is completely defined by D. It may be shown ( Backus! . Il970l ; iNorrisl . 



2006|) that 



Cijli — g \^Dij5ki + 1D k ibij — D ik 5ji — Dn5jk — Dj k 5a — Djidik 
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(A.8) 



and hence, 



(A.9) 



Finally, 



B5 = Bi + -DuDjj — 2DijDij, 



(A.10) 



from which eq. (J4]) follows. 
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